Robustness of non-Gaussian entanglement against noisy amplifier and attenuator 

environments 
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The recently developed Kraus representation for bosonic Gaussian channels is employed to study 
analytically the robustness of non-Gaussian entanglement against evolution under noisy attenuator 
and amplifier environments, and compare it with the robustness of Gaussian entanglement. Our 
results show that some non-Gaussian states with one ebit of entanglement are more robust than all 
Gaussian states, even the ones with arbitrarily large entanglement, a conclusion of direct consequence 
to the recent conjecture by Allegra et. al. [PRL, 105, 100503 (2010)]. 
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Early developments in quantum information technol- 
ogy of continuous variable (CV) systems largely concen- 
trated on Gaussian states and Gaussian operations [if . 
The sympiectic group of linear canonical transforma- 
tions [2J is available as a handy and powerful tool in this 
Gaussian scenario, leading to an elegant classification of 
permissible Gaussian processes or channels [3] . The fact 
that states in the non-Gaussian sector could offer advan- 
tage for several quantum information tasks has resulted 
more recently in considerable interest in non-Gaussian 
states, both experimental [|[ and theoretical [f|. 

Since noise is unavoidable in any actual realization 
of these information processes @ , robustness of entan- 
glement and other nonclassical effects against noise be- 
comes an important consideration. Allegra et al Q have 
thus studied the evolution of what they call photon num- 
ber entangled states (PNES) (i.e., two-mode states of the 
form \tp) = 53 c n \ n i n )) in a noisy attenuator environ- 
ment. They conjectured based on numerical evidence 
that, for a given energy, Gaussian entanglement is more 
robust than the non-Gaussian ones. Earlier Agarwal et 
al Q had shown that entanglement of the NOON state 
is more robust than Gaussian entanglement in the quan- 
tum limited amplifier environment. More recently, Nha 
et al @ have shown that nonclassical features, includ- 
ing entanglement, of several non- Gaussian states survive 
a quantum limited amplifier environment much longer 
than Gaussian entanglement. Since the conjecture of 
Ref. [?J refers to the noisy environment and the analy- 
sis in Ref. [8, 9] to the noiseless or quantum-limited case, 
the conclusions of the latter do not necessarily amount 
to refutation of the conjecture of Ref. Q. Indeed, Adesso 
has argued very recently [lfj| that the well known ex- 
tremality [11| of Gaussian states implies proof and rig- 
orous validation of the conjecture of Ref. [7| • 

In the present work we employ the recently devel- 
oped [l2[ Kraus representation of bosonic Gaussian chan- 
nels to study analytically the behaviour of non-Gaussian 



states in noisy attenuator and amplifier environments. 
Both NOON states and a simple form of PNES are con- 
sidered. Our results show conclusively that the conjec- 
ture of Ref. is too strong to be maintainable. 

Noisy attenuator environment: Under evolution 
through a noisy attenuator channel C±(k, a), k < 1, an 
input state p ln with characteristic function (CF) Xw(0 
goes to state p out with CF 



x^(6 = x^(<)^ (1 - 



2 +a)|«| 2 



(1) 



where k is the attenuation parameter [3[ . In this nota- 
tion, quantum limited channels [9| correspond to a = 0, 
and so the parameter a stands for the additional Gaus- 
sian noise. Thus, p ln is taken under the two-sided sym- 
metric action of C\ (k, a) to p out = C\ (k, a) ®C\ (k, a) (p ln ) 
with CF 



X^ t (a,6)=X^«i,< 2 )e 



-i(l-K 2 +a)(|ex| 2 + |S2| 2 ) 



(2) 



To test for separability of p out we may implement the 
partial transpose test on p out in the Fock basis or on 
Xyy te-!^)- The choice could depend on the state. 

Consider first the Gaussian case, and in par- 
ticular the two- mode squeezed state IV* (/•*)) = 
sechp X)n=o tanh™ p\n, n) with variance matrix V sq (p). 
Under the two-sided action of noisy attenuator channels 
Ci(rt,a), the output two-mode Gaussian state p out (p) — 
C\(n,a) ®C\{n,a) ( \i/>(fi)} (^(/-Ol ) nas variance matrix 



U out (^) 
V sq (ft) 



K 2 V sq (p) + (1 - 

C2 M l2 S 2/1 <73 
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where ci^ = cosh2/i, s-i^ = sinh2/i. Note that our vari- 
ance matrix differs from that of some authors by a factor 
2; in particular, the variance matrix of vacuum is the 
unit matrix in our notation. Partial transpose test [l3| 
shows that p out (p) is separable iff a > k 2 (1 — e~ 2 ''). The 
'additional noise' a required to render p out (p) separable 
is an increasing function of the squeeze (entanglement) 
parameter p and saturates at re 2 . In particular, \yj(pi)), 
pi ~ 0.5185 corresponding to one ebit of entanglement is 



2 



rendered separable when a > k 2 (1 — e^ 2 ^ 1 ). For a > k 2 , 
p° ut (p) is separable, independent of the initial squeeze 
parameter p. Thus a — n 2 is the additional noise that 
renders separable all Gaussian states. 

Behaviour of non-Gaussian entanglement may be han- 
dled directly in the Fock basis using the recently devel- 
oped Kraus representation of Gaussian channels [l2[ • In 
this basis quantum-limited attenuator Ci(/c;0), k < 1 
and quantum-limited amplifier C 2 (k;0), k > 1 are de- 
scribed, respectively, by Kraus operators [12| 



B e (n) = V m+e C e (Vl- n 2 ) e K m \m)(m + £\, 

rn—0 

i 00 , , 1 

M") = - J2 V™+ I c^(Vi-K- 2 ) e —\m + i)H: 

m— 

£ = 0, 1,2, In either case, the noisy channel 

Cj(n;a), j = 1,2 can be realized in the form £2(^2! 0) o 
Ci(/ci;0), so that the Kraus operators for the noisy case 
is simply the product set {A e i (k 2 )B£(k 1 )}. Indeed, the 
composition rule £2(^2; 0)°Ci(fci; 0) — C\{k,2K,\\ 2{k\ — 1)) 
or C2(«2Ki; 2^2(1 — k 2 )) according as K2^i < 1 or 
K2K1 > 1 implies that the noisy attenuator C±(k; a), k < 1 
is realised by the choice K2 = v/l + a/2 > 1, «i = 
k/^2 < k < 1, and the noisy amplifier £2(^50), k > 1 
by k 2 = a/k 2 + a/2 > K > 1, Ki = k/k 2 < 1 [II- -^oie 
i/iai one goes from realization o/Ci(k; a), k < 1 to that of 
£2(^5 a), K > 1 simply by replacing (l+a/2) 6y (K 2 +a/2); 
i/iis fact will be exploited later. 

Under the action of C 3 (k; a) — £2(^2! 0) oCi(«i; 0), j = 
1,2 the elementary operators |m)(n| go to 

C 2 (K 2 ;0)oC 1 (rt 1 ;0)(|m}(n|) 

00 min(m,n) 
x( K2 -V)^ + "-^(l-K 2 -W-«l)' 

x |m-£ + i)(n-^ + j|. (4) 



Substitution of k 2 = yl + a/2, «i = k/k,2 gives realiza- 
tion of Ci(k; a), k < 1 while K2 = \J n 2 + a/2, «i = k/«2 
gives that of C 2 (/c; a), « > 1. 

As our first non-Gaussian example we study the NOON 
state \ip) = (|n0) + |0n)) /v2 with density matrix 

p=i(|n)(n|®|0)(0| + |n)(0|®|0)H 

+ |0> (n\ ® \n) (0| + |0)(0| ® |n) (n| ) . (5) 

The output state p out = Ci(k;o) ® Ci(k; a){p) can be 
detailed in the Fock basis through use of Eq. (j4|) . 

To test for inseparability, we project p out onto the 
2x2 subspace spanned by the four bipartite vectors 
{|00), I On), |n, 0), \n,n)}, and test for entanglement in 
this subspace; this simple test proves sufficient for our 
purpose! The matrix elements of interest are: Pqo'oO' 



Pnn,nn, and pZ] n0 = p°^ *„. Negativity of <Ji(«,a) = 
P^OoPntnn ~ \P°0n,n \ 2 will prove for p° ut not only NPT 
entanglement, but also one-copy distillability |14 |. 

To evaluate pgg* 00 , Po„^ and P™,™, ii; suffices 
to evolve the four single-mode operators |0)(0|, |0)(n|, 
|n)(0|, and |n)(n| through the noisy attenuator Ci(n;a) 
using Eq. and then project the output to one of these 
operators. For our purpose we need only the following 
single mode matrix elements : 

xi = (n\Ci(n;a)(\n){n\)\n) 

n 

= (l + affl-^rCtftfil + a/i)-*]' 
1=0 

x [(1 - « 2 (1 + a/2)- x )(l - (1 + a/2)- 1 )]"" 1 , 
x 2 = (0|Ci(K;o)(|n)<n|)|0) 

= {l + a/2)- 1 [l- K 2 (l + a/2)- 1 ] n , 
x 3 = (0|C 1 ( K ;a)(|0)(0|)|0) = (l + a/2)- 1 , 
x 4 = (n|d( K ;a)(|0)(0|)|n) 

= (l + a/2)" 1 !!- (l + a/2)- 1 ]", 
x 5 = (n|C 1 («;a)(|n)(0|)|0) = K "(l + a/2)-(" +1 ), 

= (Old^joJdOXnDIn)* (6) 

One finds pgg* 00 = x 2 x 3 , p™* nn = X1X4, and pg£* n0 = 
x 2 /2, and therefore 

5i(k, a) = xix 2 xzXi - (\x 5 \ 2 /2) 2 , (7) 

Let cli(k) be the solution to Si(n, a) = 0. This means 
that entanglement of our NOON state survives all val- 
ues of noise a < ai(«). The curve labelled N$ in Fig.Q] 
shows, in the (a, k) space, a\{n) for the NOON state 
with n = 5: entanglement of ( 1 50} + \05})/\/2 survives 
all noisy attenuators below iVg. The straight line de- 
noted <7oo corresponds to a — k 2 : channels above this 
line break entanglement of all Gaussian states, even the 
ones with arbitrarily large entanglement. The line g\ de- 
notes a = k 2 (1 — e~ 2fil ), where fx\ = 0.5185 corresponds 
to 1 ebit of Gaussian entanglement : Gaussian entangle- 
ment < 1 ebit does not survive any of the channels above 
this line. The region R of channels above 1700 but below 
N§ are distinguished in this sense : no Gaussian en- 
tanglement survives the channels in this region, but the 
NOON state (|50) + |05))/-y/2 does. 

As a second non-Gaussian example we study the PNES 
\tp) = (| 00) + \nn)) l\pl with density matrix 

P=^(|0)(0|®|0)<0| + |0)<n|®|0)<n| 

+ \n) (0| ® \n) (0| + |n) (n\ ® |n) (n\ ) . (8) 

The output state p out = Ci(k;cl) ® Ci(k;o) (p) can be 
detailed in the Fock basis through use of Eq. (@} • 

Now to test for entanglement of p out , we project again 
p° ut onto the 2x2 subspace spanned by the vectors 
{|00), |0n), |n,0), \n,n)}, and see if it is (NPT) entan- 
gled in this subspace. Clearly, it suffices to evaluate 
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FIG. 1: Comparison of the robustness of the entanglement of 
a NOON state with that of two-mode Gaussian states under 
the two-sided action of symmetric noisy attenuator. 



FIG. 3: Comparison of the robustness of the entanglement 
of a NOON state with that of all two-mode Gaussian states 
under the action of two-sided symmetric noisy amplifier. 
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FIG. 2: Comparison of the robustness of the entanglement of 
a PNES state with that of two-mode Gaussian states under 
the action of two-sided symmetric noisy attenuator. 



the matrix elements pg^ „, p°o*„ , and pfj$ nn , for if 
6 2 ( K ,a) = /C* „/Co!„o - \p° £ nn \ 2 is negative then p° ut 
is NPT entangled, and one-copy distillable. 

Once again, the matrix elements listed in prove suf- 
ficient to determine <5 2 (k, a): Pg"!n0 = PnanO = ( x i x 2 + 
x 3 x 4 )/2, and p^ nn = \x 5 \ 2 /2, and so 

8 2 {K,a) = {(x lX2 + x 3 x 4 )/2) 2 - (\x 5 \ 2 /2) 2 . (9) 

Let d2( K ) denote the solution to (^(k, a) = 0. That 
is, entanglement of our PNES survives all a < 02 (k). 
This ci 2 (k) is shown as the curve labelled P5 in Fig. [2] for 
the PNES (| 00) + \55))/V2. The lines g Y and 5oo have 
the same meaning as in Fig.[T] The region R above g^ 
but below P5 corresponds to channels (k, a) under whose 
action all two-mode Gaussian states are rendered sep- 
arable, while entanglement of the non-Gaussian PNES 
(| 00) + |55))/V2 definitely survives. 

Noisy amplifier environment: We turn our attention 
now to the amplifier environment. Under the symmetric 
two-sided action of a noisy amplifier channel C2 (k; a) , k > 
1, the two-mode CF Xw(£i>&) 1S taken to 

x ^(Ci,6) = x^«i,< 2 ) e ^ (K2 - 1+a)(l « l|2+l « 2|2) . 



In particular, the two-mode squeezed vacuum state 
\ip[fj,)} with variance matrix V sq (/i) is taken to a Gaussian 
state with variance matrix 

V out (p) = n 2 V sq {p) + (k 2 -1 + a)l 4 . (f0) 

The partial transpose test [l3[ readily shows that the out- 
put state is separable when a > 2 — k 2 (1 + e~ 2/J ): the 
additional noise a required to render the output Gaus- 
sian state separable increases with the squeeze or en- 
tanglement parameter \x and saturates at a = 2 — k . 
for a > 2 — k 2 the output state is separable for ev- 
ery Gaussian input. The noise required to render the 
two-mode squeezed state \ip(pi)} with 1 ebit of entangle- 
ment (fix sa 0.5185) separable is a = 2 — n 2 (l + e~ 2pi ). 
Now we examine the behaviour of the NOON state 
(|n0) + \0n))/y/2 under the symmetric action of noisy 
amplifiers C2(^;a), k > 1. Proceeding exactly as in the 
attenuator case, we know that p out is definitely entan- 
gled if 5 3 {K,a) ee pg^bo/CVn - |Po"!nol 2 is negative. As 
remarked earlier the expressions for C\(n;a), k < 1 in 
Eq. ([6]) are valid for C 2 (/t;a), k > 1 provided 1 + a/2 is 
replaced by n 2 + a/2. For clarity we denote by 2; ■ the ex- 
pressions resulting from Xj when Ci(k; a), k < 1 replaced 
by C 2 (k; a), > 1 and 1 + a/2 by k 2 + a/2. For instance, 
x' 5 ee (n\C 2 (n;a)(\n)(0\)\0) = n n (n 2 + a/2)- {n+1 ^ and 
5 3 (k; a) = x' lX2 x'zx' 4 - (|x5| 2 /2) 2 . 

Let 0,3(0) be the solution to Ss(k, a) — 0. This is rep- 
resented in Fig. [3] by the curve marked N5, for the case 
of NOON state (|05) + |50))/V2. This curve is to be 
compared with the line a = 2 — k 2 , denoted g^, above 
which no Gaussian entanglement survives, and with the 
line a — 2 — K 2 (l + e~ 2Al1 ), pi = 0.5185, denoted g\, above 
which no Gaussian entanglement < 1 ebit survives. In 
particular, the region R between g^ and N$ corresponds 
to noisy amplifier channels against which entanglement 
of the NOON state (|05) + |50))/V2 is robust, whereas no 
Gaussian entanglement survives. 

Finally, we consider the behaviour of the PNES ( 1 00) + 
\nn))/\/2 in this noisy amplifier environment. The out- 
put, denoted p out , is certainly entangled if 64 (k, a) = 
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FIG. 4: Comparison of the robustness of the entanglement 
of a PNES state with that of all two-mode Gaussian states 
under the action of two-sided symmetric noisy amplifier. 

Pon!onPnO%o " IPoofmJ 2 is negative. Proceeding as in the 
case of the attenuator, and remembering the connec- 
tion between Xj J s and the corresponding Xj's, we have 

<5 4 (k, a) = + x^x' 4 )/2) 2 - (\x' 5 \ 2 /2) 2 . The curve 

denoted P5 in Fig. 2] represents 04 (re) forming solution to 
<5 4 (re,a) = 0, for the case of the PNES (|00) + |55))/V2. 
The lines and g\ have the same meaning as in Fig. [3] 
The region R between g^ and P5 signifies the robust- 
ness of our PNES : for every re > 1, the PNES is seen to 



endure more noise than Gaussian states with arbitrarily 
large entanglement. 

We conclude with a pair of remarks inspired by in- 
sightful comments by the Referee. First, our conclu- 
sion following Eq. ^ and Eq. (|10l) that entanglement 
of two-mode squeezed (pure) state \ip (//)) does not sur- 
vive, for any value of /i, channels (re, a) which satisfy 
the inequality |1 — k 2 \ + a > 1 applies to all Gaussian 
states. Indeed, for an arbitrary (pure or mixed) two- 
mode Gaussian state with variance matrix Vq it is clear 
from Eqs. ([3)), (fTO]) that the output Gaussian state has 
variance matrix V out = k 2 V g + (|1 - k 2 \ + a)l 4 . Thus 
|1 — re 2 | + a > 1 immediately implies, in view of nonnega- 
tivity of Vq, that V out > I4, demonstrating separability 
of the output state for arbitrary Gaussian input [l3| . 

Secondly, Gaussian entanglement resides entirely 'in' 
the variance matrix, and hence disappears when environ- 
mental noise raises the variance matrix above the vacuum 
or quantum noise limit. That our chosen states survive 
these environments shows that their entanglement resides 
in the higher moments, in turn demonstrating that their 
entanglement is genuine non-Gaussian. Indeed, the vari- 
ance matrix of our PNES and NOON states for N = 5 is 
six times that of the vacuum state. 
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